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. We derive an evolution equation describing the high energy behavior of the cross section 

for the single diffractive dissociation in deep inelastic scattering on a hadron or a nucleus. 
The evolution equation resums multiple BFKL pomeron exchanges contributing to the cross 
Q\ ■ section of the events with large rapidity gaps. Analyzing the properties of an approximate 

solution of the proposed equation we point out that at very high energies there is a possibility 
"*^h , that for a fixed center of mass energy the cross section will reach a local maximum at a certain 

intermediate size of the rapidity gap, or, equivalently, at some non-zero value of the invariant 
I ' mass of the produced particles. 

I. INTRODUCTION 

Some time ago an evolution equation was derived in the framework of Mueller's dipole model [|l]~[5} which 
resums multiple Balitsky, Fadin, Kuraev and Lipatov (BFKL) pomeron exchanges for deep inelastic 
scattering (DIS) in the leading In 1/x approximation 0. In terms of the usual Feynman diagram formulation 
the equation in sums up the so-called pomeron "fan" diagrams contributing to the total cross section 
of a quark-antiquark pair on a hadron or nucleus (see Fig. [j]), and turns out to be a generalization of the 
Gribov, Levin and Ryskin (GLR) equation Similar nonlinear equation for multiple pomeron exchanges 
was also obtained by Balitsky || in the effective lagrangian approach. 

The physical picture of DIS presented in |7j is the following: virtual photon splits into a quark-antiquark 
pair, which, by the time it reaches the target develops a cascade of dipoles, each of which independently 
interacts with the target. All the QCD evolution was included in the developed cascade of the dipoles in the 
large N c limit, similar to 

The equation in Jt| was written for the elastic amplitude of the scattering of the qq pair of transverse 
size x located at impact parameter b on a hadron or nucleus, N (x, b, Y), where Y ~ In 1/x is the rapidity 
variable of the slowest particle in the qq pair. The target's structure function F2, as well as the total DIS 
cross section, could be easily expressed in terms of -/Vq(x, b, Y): 
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where the quark and antiquark are located at transverse coordinates Xo and Xi correspondingly and xoi = 
xi — xq. Q 2 is the virtuality of the photon, z is the fraction of the pair's momentum carried by a quark. 
$(xoi,z) is the probability of virtual photon fluctuating into a qq pair, which we will refer to as virtual 
photon's wave function. The exact expression for <!>(xoi, z) is well known and is given in [Q as well as in a 
number of other references. 

In the following equation was written for iV (x, b, Y): 
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where the initial condition 7(xoi,b) was given by the Glauber-Mueller interaction formula of the quark- 
antiquark pair with the nucleus 



7(x i,b ) = 1 - exp 



2N C S ± 



(3) 



Here, for a cylindrical nucleus, S± = ttR 2 is the transverse area of the hadron or nucleus, A is the atomic 
number of the nucleus, and xG is the gluon distribution in a nucleon in the nucleus, which was taken at 
the two gluon (lowest in a) level [ fL0|| . In Eq. (||) p is an ultraviolet regulator, which never appears in 
JV (xoi,b,y). 




FIG. 1. An example of the pomeron "fan" diagram, which were considered in deriving (g). 

There have been several attempts to solve Eq. (^|) analytically Jl^Jl^]. Unfortunately it turned out to be 
very hard to construct an analytical solution describing simultaneously the behavior of A/o(xoi,b, Y) both 
outside and inside of the saturation region. Instead Eq. (Q) was solved separately outside [k > Q S (Y)] and 
inside [k < Q S (Y)] the saturation region (l^Jl^], where approximately |ll[] 



Q a {Y) ~ A a 2 A 1 / 3 
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with dp - 1 = iaC F j n 2 j n j[l],[L4| it was concluded that for not very large energies, corresponding to 
rapidities of the order of Y ~ 1/a (or, equivalently, for k > Q S {Y) ), the solution of Eq. (||) behaves with 
energy like a single BFKL pomeron exchange, with multiple pomeron corrections setting in as energy increases 
and slowing down the growth of 7Vo(xoi, b, Y). For very high energies (Y > a ^_ 1 In or k < Q S (Y)) the 
solution of Eq. (|2|) saturates to a constant, namely iVo(xoi, b, Y) = 1, which corresponds to the blackness 
of the total cross section of the quark-ant iquark pair on the nucleus Thus, as one can see at least 

the qualitative features of N (x±,h,Y) as a function of Y and x± arc known. For a better quantitative 
understanding of the behavior of JVrj(xj_, b, V) one should, probably, solve Eq. (||) numerically, which would 
be a very interesting and useful project to perform. 




FIG. 2. Single difTractive dissociation process considered in the paper. The interaction between the target and 
virtual photon is represented schematically by an exchange of a color singlet object. 

In JtJ the object of interest was the total inclusive DIS cross section, where no restrictions are imposed on 
the final state of the process. In this paper we are going to study the cross section of the single diffractive 
dissociation. The physical picture of the process we are going to consider is the following: in DIS the virtual 
photon interacts with the hadron or nucleus breaking up into hadrons and jets in the final state. At the same 
time the target hadron (nucleus) remains intact. The particles produced as a result of hadron's breakup do 
not fill the whole rapidity interval, leaving a rapidity gap between the target and the "slowest" produced 
particle. The process is depicted in Fig. g. In diffractive dissociation one imposes a restriction on the final 
state — the existence of a rapidity gap. 




FIG. 3. Traditional description of single diffractive dissociation. Dash-dotted line represents the final state. 

Below we are going to employ the techniques of Mueller's dipole model similarly to the way they 

were applied in M, to construct a cross section of the single diffractive dissociation in DIS which would 



3 



include the effects of multiple pomeron exchanges. Analogous to [0,0 we will neglect the diagrams with 
pomeron loops, i.e. , the diagrams where a pomeron splits into two pomerons and then the two pomerons 
merge into one pomeron. As was argued in [l3| these diagrams are suppressed compared to the pomeron 
fan diagrams of Fig. |] in DIS on a large nucleus by factors of A 1 / 3 . Thus one can safely neglect them up 



to rapidities of the order of Y > 



■ In ■ 



For the case of DIS on a proton, strictly speaking there is 



ap — 1 a 2 A 1 / 3 ' 

no similar argument allowing one to neglect pomeron loop diagrams. Therefore, whereas the fan diagram 
evolution dominates in DIS on a nucleus, it could only be considered as a model for DIS on a hadron. 

In the traditional description of the single diffractive dissociation one usually considers triple pomeron 
vertex [p|JT^ where the pomeron above the vertex is cut, and the two pomerons below the vertex are 
on different sides of the cut, thus producing a rapidity gap, as shown in Fig. ^. That way the particles are 
produced by the cut pomeron in the rapidity interval adjacent to the virtual photon's fragmentation region 
and no particles are produced with rapidities close to the final state of the target due to uncut pomerons. 
In this paper we want to enhance this picture by including multiple pomeron exchange diagrams. We would 
like to understand the behavior of the diffractive dissociation at very high energies, which may include some 
effects of saturation of hadronic or nuclear structure functions. 

An example of a graph in the usual Feynman diagram language which we will consider below is given in 
Fig. ||. Dash-dotted line corresponds to the final state, i.e., to the cut. As one can see in Fig. |some of the 
pomerons are cut, some remain uncut. In the region of rapidity where pomerons are cut we have particles 
being produced. In the notation of Fig. ^ this corresponds to the interval in rapidity from Y — Y a to Y. 
In the region where the pomerons are uncut (rapidity interval from to Yq) nothing is produced, which 
corresponds to a rapidity gap. That way Fig. demonstrates a generalization of the traditional picture of 
Fig. (H), in which the cut pomeron splits into two pomerons, which later branch into two uncut pomerons 
each. Our goal in this paper is to resum all the diagrams where the cut pomeron can split into any number 
of cut pomerons via fan diagrams, and the cut pomerons in turn split into uncut pomerons, which can also 
branch into any number of uncut pomerons interacting with the target below. 





FIG. 4. A diagram contributing to the diffractive dissociation with rapidity gap Yo as considered in the text. 



In Sect. II we employ the dipole model formalism [[!]-§} to write down an equation for the cross section 
of the single diffractive dissociation of a quark-antiquark pair on a hadron or nucleus with the rapidity gap 
bigger or equal to Yo- The task is not very straightforward, since the dipole model provides us with the 
dipole wave function of the qq pair by the time it hits the target. In the rest frame of the target the typical 
time scale associated with the interaction is negligibly small compared to the typical lifetime of the QCD 
evolved wavefunction of the pair. Thus one could say that the interaction is instantaneous and happens at 
the time t — |Q. However, after the interaction, the qq state with many gluons in it can undergo several 
transformation until all the partons in it reach the mass shell at t — oo. Dipole model does not provide us 
with any information about these final state interactions. Nevertheless one can gain control over the final 
state through the unitarity condition and cancelation of certain classes of final state interactions, as will be 
shown below. As a result we obtain a nonlinear evolution equation for the diffractive cross section, shown 
in formula (M). 
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We show how Eq. (|9|) can be re-derived by applying the AGK cutting rules to the pomeron fan diagrams 
in Sect. III. Since the cut or uncut pomerons precisely specify the final state in the traditional language, we, 
therefore, proved the consistency of our treatment of the final states in the dipole model. 

Even though the exact analytical solution of that equation seems to be very complicated in Sect. IV we 
analyze the properties of the solution. As a toy model we consider a simplified version of the equation with 
the kernel independent of transverse coordinates. We show that the simplified equation's solution exhibits 
a remarkable property at very high energies. We plot the cross section of diffractive dissociation events 
with rapidity gap Yq at fixed large center of mass energy as a function of the size of the rapidity gap Yq . 
This is equivalent to plotting the cross section as a function of the invariant mass of the produced particles 
M|, since, as could be easily seen from Fig. ^, Yq = Y — lnMj^/Q 2 ''. At not very large Yq (large M^) 
the cross section increases with the increase of rapidity gap, which agrees with the result of the traditional 
triple pomeron description of the diffractive dissociation of pjl5[]. However, as Yq gets very high comparable 
to the values of rapidity at saturation, which corresponds to a very small produced mass M|, the cross 
section reaches a maximum and starts decreasing. That way if one would be able to measure the single 
diffractive cross section for a range of different invariant masses of the produced particles, one should be able 
to observe the maximum and turnover of the cross section for small masses if the energy is high enough for the 
saturation effects to become important. If the exact, probably numerical, solution of the equation confirms 
our conclusion, which was reached by approximate methods, then the single diffractive cross section is a new 
and independent observable which can be used to test whether the saturation region has been reached in 
DIS experiments. 

The summary of our results is presented in Sect. V. 



II. EVOLUTION EQUATION FOR DIFFRACTIVE CROSS SECTION 



In this section we will write down an evolution equation resumming multiple pomeron exchanges for the 
cross section of single diffractive dissociation in DIS on a hadron or nucleus (see Figs. || and ||). As was 
mentioned in the introduction the task is complicated by the fact that the techniques of the dipole model 
which we are going to employ provide us with the state of the system at time t = 0, without putting 
any restrictions on its subsequent time evolution until t = oo. The dipole cascade is developed by the time of 
interaction and interacts with the nucleus at t = for any process. It is the subsequent time evolution of the 
cascade that determines whether the event is going to have a rapidity gap, which would correspond to some 
of the dipoles recombining, or not. This evolution is not included in the dipole model However we 

will be able to constrain that time evolution after the interaction took place and to construct the diffractive 
cross section using the following two observations. 




FIG. 5. An example of elastic scattering of the quark-antiquark pair on the nucleus as pictured in the text. 



The first observation is that if one knows the total cross section for a the scattering of some probe on 
a target one can easily construct the elastic scattering cross section out of it using unitarity constraints. 
Really, if the total cross section is |19 



atot = 2 / d 2 b [1 - S(b)} 



(5a) 
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with S(b) the S'-matrix of the scattering process at a fixed impact parameter b, then the elastic cross section 
is 



d 2 b [1 



S(b)f 



(5b) 



Using Eqs. (EaJ),(5b) it was shown in |fl 19 that the elastic cross section for the scattering of the quark- 
antiquark pair generated by a virtual photon in DIS on a nucleus or hadron could be expressed in terms of 
the elastic amplitude iVo (xoi , b, Y), or, equivalently, in terms of the imaginary part of the forward scattering 
amplitude iVo(xoi, b, Y)d 2 b (The forward amplitude is purely imaginary). The result was [ p~T]JTg| l 



iV e '(x i,b,F) = 7V (x 01 ,b,y) 2 



(0) 



and is illustrated in Fig. g. The cut in Fig. [| corresponds to the final state, when the particles reached 
the mass shell and t = oo. Thus in that picture we explicitly impose the constraint that not only the qq 
system developed a dipole cascade and the cascade in turn interacted with the target at t = 0, but it also 
recombined back into a qq pair by the time the system reached the final state at t = oo. Eq. (g) allows us 
to calculate the elastic cross section of the scattering of any (not necessarily original) dipole in the dipole 
cascade on the nucleus including multiple pomeron evolution. 

The second observation which we have to make before starting to derive our equation is based on the work 
of Chen and Mueller H and concerns the so-called final state interactions. By that we mean the interactions 
such as branching or recombination of partons in the dipole cascade after t = 0, i.e., after the interaction 
with the target. Chen and Mueller Q succeeded in proving that certain classes of final state interactions 
cancel (see Fig. ||). 








B 







FIG. 6. Classes of diagrams with final state interactions which cancel. 

We denote a gluon by the double line in Fig. || corresponding to the large N c limit. In each diagram of Fig. 
H we assume that the summation over all possible connections of the gluon to the quark and antiquark lines 
is performed. In Fig. we consider the situation when the gluon is emitted in a dipole after the interaction 
with the target at t = 0, which is denoted by dotted line. The three different positions of the t = oo final 
state cut, denoted by dash-dotted line in Fig. ||A correspond to the cases when the emitted gluon is present 
in the final state and when the gluon recombines back in the amplitude or complex conjugate amplitude 
leaving the dipole intact by t ~ oo. As was proved in the sum of the three cuts of Fig. is zero. 

In Fig. |^B we consider the situation when the gluon was developed in the dipole wave function before 
the interaction with the target and is already there by the time t = 0. Then it can either remain in the 
final state t = oo or it can recombine back, leaving only the original dipole in which it was produced in the 
final state. In Fig. |^B we explore the case when the gluon in the final state becomes a gluon emitted after 
t = in the complex conjugate amplitude (first diagram in Fig. ^3). Chen and Mueller Q showed that this 
diagram is canceled by the second diagram in Fig. ^|B. 

Diagrams of Fig. |6| demonstrate that the contribution of gluons produced or absorbed after t = cancels 
out in calculation of inclusive quantities such as the total cross section. We can now conclude that the total 
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inclusive DIS cross section can be calculated using just t = formalism, as was done in |7]]. Really, as one 
can see from Fig. ^| all the final state interactions cancel, effectively leaving the t — state unchanged by 
t = oo. This conclusion is intuitively easy to understand: one just needs to have some interaction of the 
projectile with the target to include the diagram in the total cross section, independent of what happened 
after that interaction. 

Now we are ready to write down our equation for diffractive cross section. Let us start by defining the 
object for which the equation will be written. We denote by N D (x±, b, Y, Yq) the cross section of single 
diffractive dissociation of a dipole of transverse size xj_, rapidity Y and impact parameter b on a target 
hadron or nucleus. The process has a rapidity gap covering a rapidity interval adjacent to the target (see Fig. 
|,|] ) which could be greater or equal than Yq. The corresponding diffractive structure function, similarly to 
Eq. is 

Fi D ( x ,Qi,Y ) = -^- I f^^ X01)Z ) d H 7v D (x 01 ,b,y,y ). (7) 

4:7T OiEM J 47T 

If one wishes to obtain a cross section of diffractive dissociation of the dipole with a fixed rapidity gap Yq 
one has to differentiate N D with respect to Yq, as will be discussed later. When Y = Yq the rapidity gap 
fills out the whole rapidity interval, turning theprocess of dipole's dissociation into an elastic scattering. 
Taking the answer for elastic process from Eq. (|6|) we obtain an initial condition for the evolution of N D 

iV D (x ± ,b,y =Y ,Y ) =iV 2 (x ± ,b,r ). (8) 

Note that Eq. (||) insures that after the interaction with the target at t = the system of developed 
dipoles only recombines back into the original dipole by t — oo. All other possible final state fluctuations 
have been excluded because we explicitly put Nq (xj_, b, Yq) as the amplitude of the elastic process. 

To construct N D one has to require that nothing is produced in the t = oo final state in the rapidity 
interval from to Yq, where the target hadron or nucleus is situated at rapidity 0. This does not restrain the 
rapidity gaps from being greater than Yq. That way N D would include all events with rapidity gap greater 
or equal to Yq. The condition of the rapidity gap up to Yq is satisfied by Eq. (||). 

As was mentioned above our physical picture of the diffractive dissociation event is the following: before 
hitting the target the virtual photon develops a dipole cascade, which at the time t = interacts with the 
target. After that some partons in the cascade may recombine, producing rapidity gaps, some may not. By 
imposing the recombination condition in Eq. (|8|) we made sure that the dipoles with rapidities y < Yq will 
recombine by t — oo, producing a rapidity gap. The dipoles with y > Yq are free to either recombine into 
the dipoles off which they were produced or to remain unchanged till t — oo. As long as y > Yq there are 
no restrictions on the final state like rapidity gaps. Therefore the dipoles are free to recombine after t = 0. 
However, for that situation the cancelation rules proven by Chen and Mueller Q apply. The graphs when 
the dipoles recombine to either increase the size of the rapidity gap Yq or to produce extra rapidity gaps 
at y > Yq are included, but as one can see from Fig. they cancel. For the asymmetric case when the 
dipole in consideration interacts with the target in the amplitude and then recombines back into the dipole 
off which it was produced which has no interaction with the target in the complex conjugate amplitude we 
employ a different cancelation mechanism: the contributions when the dipole recombines in the complex 
conjugate amplitude at the times t > and t < cancel each other, leaving us with the diagram of similar 
to the second graph of Fig. ^|B, which could be incorporated in the elastic interaction of the original dipole 
with the target A^o(xj_, b, y). That way the asymmetric graphs may include the dipoles of any rapidities 
y > Yq interacting elastically with the target via Eq. (ph. For the "symmetric" case when the dipole at 
y > Yq interacts elastically with the target in the amplitude and then recombines into the original dipole 
which interacts with the target elastically in the complex conjugate amplitude we applied the cancelation 
of Fig. ^|B. Thus the symmetric case reduces to the situation when the same dipole interacts elastically in 
the amplitude and complex conjugate amplitude at y — Yq and is put only in the initial condition of our 
equation given by Nq (xj_, b, Yq). 

Dipole splitting is a little more complicated. The softest dipoles in the cascade interact the way it is 
described in Eq. (^), which fixes the final state for them. Now let us consider a dipole which is not one of 
these soft dipoles and, therefore, does not interact with the target. The dipole can split into two dipoles 
either before or after t = 0. This event is allowed only when the produced real gluon has the rapidity greater 
or equal than Yq, so that it would not destroy the rapidity gap. In that case the real gluon emitted before 
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t = is canceled by a virtual correction in the same dipole before t = 0. The real gluon emitted after 
t = is canceled by the virtual term after t = 0, similarly to Fig. ||A. Finally, this non-interacting dipole 
can develop a virtual corrections with the gluon being softer than Yq, i.e. having smaller rapidity than the 
rapidity of the gap. These gluons could not be canceled by the real contributions, since those are excluded 
by the condition of the gap's existence. However, since the virtual corrections could happen either before or 
after t = 0, these two combinations come with different signs. The absolute values of both t > and t < 
contributions in that case are the same and, therefore, they cancel each other. Thus we have shown that the 
final state dipole splitting is canceled together with the virtual corrections in the non- interacting dipoles. 

The above discussion is summarized in Fig. 0. The change of the cross section of diffractive dissociation 
N D in one step of the rapidity evolution could be generated by several terms on the right hand side of the 
equation in Fig. [t]. The first terms corresponds to the usual single BFKL pomeron evolution, while the rest 
of the terms correspond to the triple pomeron vertex, similarly to the equation written in |l|,|7| . Effectively 
Fig. [t] states that in one step of evolution the color dipole splits into two dipoles. The subsequent evolution 
can happen in one of the dipoles, which corresponds to the first term on the right hand side of Fig. [?]. It 
could also happen that we continue the evolution leading to the cross section of diffractive dissociation in 
both dipoles (the second terms in Fig. [?]). The third and forth terms correspond to the asymmetric cases. 
As was discussed above in the asymmetric case the dipole can interact elastically with the target in the 
amplitude but not in the complex conjugate amplitude and vice versa at any rapidity Y > Yq. Thus we 
include the cases when only one of the produced dipoles interacts asymmetrically with the target (the third 
term in Fig. 0) and when both of them do so (the fourth term in Fig. |7|) . The coefficients in front of different 
terms come from combinatorics. Note that they match the coefficients of the Abramovskii, Gribov, Kancheli 
(AGK) |L6| cutting rules for two pomeron contributions. 





FIG. 7. The equation for the cross section of diffractive dissociation. 

Recalling the initial condition of Eq. (@) we can turn the picture of Fig. [7] into a formula. For that we also 
have to include the virtual corrections similarly to the way they were included in jlj . The resulting equation 
is 

WV b, Y Y ) = iV 2 (x 01 , b, Yo) h ffl^o) + Tdy KT)^) 

x f d 2 x 2 I 2 N ° ( x 02 , b + ic 12 , y , Y ) + N D (x 02 ,h + ^ Xl2 ,y,Y ) N D (x 12 , b + ix 02 , y, Y ) 



-4A^ D (x 02 ,b+ -x 12 ,y,yo)^o(xi 2 ,b+ -x 02 ,y) + 2A (x 02 ,b+ -x 12 ,y) A^ (x 12 ,b+ -x 02 ,y)]. (9) 

The shifts in the impact parameter dependence of the functions on the right hand side of Eq. (^) have 
occurred because the centers of mass of the produced dipoles are located at different impact parameters 
than the center of mass of the initial dipole. 

Eq. ([)]) describes the small-x evolution of the cross section of the single diffractive dissociation for DIS 
N and is the central result of this paper. 
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III. REDERIVING EVOLUTION EQUATION WITH AGK CUTTING RULES 



Eq. <Mj can also be derived using the technique employing the AGK cutting rules . In this section we 
will not employ the dipole model. Instead we will consider pomeron fan diagrams in the usual perturbation 
theory, imposing a restriction that there are only triple pomeron vertices in the theory. This assumption is 
true only in the large N c limit, as was shown in the dipole model. 

AGK cutting rules jl6| apply to the BFKL pomerons in QCD (sec |^2| for a detailed discussion). For the 
BFKL pomeron one can write 



-.BFKL 



is the one BFKL pomeron contribution to the forward scattering amplitude and G 



(10) 



BFKL 



IS 



where 

the total inelastic cross section. In writing Eq. ( |Toj ) we have neglected the elastic contribution to the cross 
section, which is a justified assumption for the BFKL pomeron exchange f2i| . Eq. (10) tells us that the 
inelastic cross section in one pomeron exchange approximation is twice the forward amplitude of the process, 
since the latter is purely imaginary. 

ABC 





4N D N 




FIG. 8. Different pomeron cuts contributing to the cross section of diffractive dissociation which lead to different 
terms on the right hand side in Eq. (bl). 



Different cuts of pomeron fan diagrams corresponding to the terms on the right hand side of Eq. (^) are 
displayed in Fig. M. The first (linear) term on the right hand side of Eq. (m corresponds to just usual 
linear BFKL evolution with the factor of 2 arising from Eq. (|l^) and is not shown in Fig. pi The non- 
linear evolution corresponds to a pomeron splitting into two via the triple pomeron vertices. Let us now 
concentrate on the first (topmost) triple pomeron vertex in the diagrams of Fig. ^ This is where one step 
of non-linear evolution takes place in our picture. The graph of Fig. ||A corresponds to the case when a cut 
pomeron splits into two cut pomerons, which later on produce rapidity gaps. That is the cross section of 
diffractive dissociation Nd splits into two similar cross sections. Therefore this term gives rise to N% ter m 
in Eq. (^|) . The factor in front of this term results from multiplying 2 from the cut pomeron on top of the 
vertex by the symmetry factor of 1/2. The diagram of Fig. |^B corresponds to a cut pomeron splitting into 
another cut pomeron, which then develops a rapidity gap (Njy), and an uncut pomeron, corresponding to 
the forward scattering amplitude No- Since there are two of such diagrams, multiplying the factor of 2 of the 
cut pomeron by 2 and putting a minus sign since one of the resulting pomerons is completely on one side of 
the cut we will get —4Ne>Nq, which is the contribution of that diagram in Eq. (|). Finally the graph of Fig. 
||C corresponds to the case when both of the pomerons below the vertex are uncut. Each uncut pomeron 
gives a factor of Nq. The factor of 2 arises from the cut pomeron above [Eq. (|Tc|)j_. 

Therefore we were able to reproduce all the terms on the right hand side of Eq. (m from the usual pomeron 
fan diagrams using the AGK cutting rules. We note that a similar result was obtained in the framework of 
reggeon theory in p3|. The cut (or uncut) pomeron specifies t = oo final state very precisely: if the pomeron 
is cut it insures that there will be particles produced in the final state. Uncut pomeron insures that nothing 
is produced in the t — oo final state. Thus in this section we have demonstrated the consistency of our 
analysis of the evolution of the dipole state after the interaction with the target which was employed in the 
previous section. 
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IV. PROPERTIES OF THE SOLUTION OF THE EVOLUTION EQUATION 



Eq. (^|) can be rewritten in the differential form. Similar to what was done in M we differentiate both 
sides of Eq. with respect to Y and once again make use of Eq. (^) to obtain 



ajy^xoi.b, y,y ) 2 a c F . 2 
qy = — / d ** 



1 01 



2 2 
X 02 X 12 



2nS (xqi - X02) In 



A^(x 02 ,b+~x 12 ,y,Y ) 



a ° F f d 2 x 2 ^ 



f\ [iv D (x 02! b + i Xl2 , y y ) iv D (x 12 , b + lx 02 , y y ) 

• Z '02 :1 '12 " * 



4 7V D (x 02 , b + ixia, y y ) JV (xi2, b + ix 02 , y) + 2iV (x 02 , b + lxi2, y) ^o(xi2, b + 1 X0 2, Y)], (11) 



with Eq. (g) being the initial condition for the differential equation (|ll| 
Let us define the following object 



^(xoi, b, y y ) = 2 jvo(xoi, b, y) - a^(x i, b, y y ) 0(y - y ), 



(12) 



which has the meaning of the cross section of the events with rapidity gaps less than Yq. The theta-function 
that multiplies N D in Eq. ([l2|) insures the trivial fact that the rapidity gap can not be larger than the total 
rapidity interval. One can see that Eq. (O) can be rewritten for Yq < Y in terms of T as 



&F(xoi,b,Y,r ) _ 2aCV / , 2 



dY 



d x 2 



X 02 X 12 



27rJ 2 (x i - x 02 ) In 



J r (x 02 ,b + lxi2, yy ) 



aC F 



d 2 x 2 



'-01 



1 



1 



J^(x 2,b+ -xi2,y y ) J r (xi2,b+ ~x 02 ,Y,Yo), 



02-^12 



with the initial condition 

•F(xoi, b, y = y , y ) = 2iv (x i,b, y ) - jv 2 (xoi, b, y ). 



(13) 



(14) 



If one differentiates Eq. (||) with respect to Y one would obtain the same differential equation as Eq. ( |i"3| ) 
for N 



0JV o (xoi,b,Y) 2aC F f , 2 
= ^T- M»a 



XnoXi 



27rJ (xqi - X02) In 



^01 



iV (x 2,b+ -xi 2 ,y) 



aCp 



d 2 x 2 



X 2 1 1 

-g-22g_ A (x 02 , b + -X12, Y) iV (xi2, b + -x 02 , Y). 

X 02 X 12 ^ * 



(15) 



(16) 



The initial condition for Eq. (|T^) is different from that for Eq. ( |l~5| ) 

JV (xoi,b,y = 0)= 7 (x i,b) 

with 7(xoi, b) given by Eq. (||). From Eq. ( |l2| ) one can see now that in order to find iV D (xoi, b, Y, Yq) one 
has to solve the same equation [(|l3|) or Jl5|)] for two different initial conditions given by formulae (|l4|) and 
(|l6|), which would yield us with the results for N and T . Then N D could be found using Eq. (12). 

As was mentioned above the exact analytical solution of Eq. ([l3|), or, equivalently, (]l5|) has not been 
found jfl 14|. Nevertheless the qualitative behavior of the solution of this equation is understood very well, 
with some quantitative estimates performed |Ll],[l4|. It turned out that the solution of this equation can 
be qualitatively well approximated by the solution of the equation with suppressed transverse coordinate 
dependence 
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=(a P -l)N (Y)-(ap-l)N (Y) 2 . (17) 



dY 



Eq. (|l7|) is a toy model of the full Eq. ( |l5| ) , with all the transverse coordinate integrations suppressed and 
the dipole BFKL kernel substituted by its eigenvalue at the BFKL saddle point. The fact that the coefficients 
in front of the linear and quadratic terms in Eq. ( |j~7| ) are identical could, for instance, be understood in the 
double logarithmic limit (^02,^12 <§C Xqi), in which the kernel in front of the quadratic term in Eq. ( |l5| ) 
produces an extra factor of two and becomes equal to the kernel in front of the linear term (l^Jl^l ■ However 
it is a more general property of the solution of Eq. (|l5|) , which is valid not only in the double logarithmic 
limit. 

In this paper we will try to understand the qualitative behavior of the solution of Eq. ([)]) using the toy 
model of Eq. (jl7|). An exact, probably numerical, solution of Eq. ( |T5| ) would give more precise results 
for N D . However, we believe that the qualitative features provided by the toy model of Eq. ( )l7| ) will be 
preserved in the exact solution. 

Let us suppose that Nq(Y) is a solution of Eq. (fjl|), or, in our toy model, of Eq. (p"7|). Then T{Y, Yq) will 
be a solution of the same equation with different initial condition set at a different value of rapidity (see Eq. 
(jjlj)). Thus, if we neglect the transverse coordinate dependence, the solution for T{Y, Yq) could be obtained 
from the solution for Nq(Y) just by a shift in rapidity. Define the shift variable AY" by 

No(Yo+AY) = 2Nq(Yq)-Nq(Yq) 2 (18) 

so that 

T(Y,Y Q ) =N Q (Y + AY). (19) 



Now one can see that since No(Y) satisfies the differential equation ( |15| ) (or Eq. ([17]) in our toy model), 
then T{Y,Yq) given by Eq. (19) satisfies the same equation ( jl7| ) with the initial condition of Eq. (|l4|). Of 



course this is true only when one neglects the transverse coordinate dependence of N and T, but could also 
be a reasonable approximation for the case of weak transverse coordinate dependence. 
Solving Eq. ( p"7j ) with the initial condition No(Y = 0) = 7 we obtain 

Nq{Y) = '—, TW r, (20) 

uv ' 1 + 7 (e( Q ^- 1 ) y - 1) v 1 

which employing Eq. (jlj) yields us with the following expression for AY" 

AY = — In (2 + — — e ( ap - 1)Yo ). (21) 



ap — 1 \ 1 — 7 
Eq. ( |l9| ) allows us to write down a toy model expression for T 

Hy, Y ) = ^~r^ ( 22 ) 

1 _ 7 + 7e (a P -i)y (2 + _L_ e { ap -i)Y a \ 
Now we can explore the qualitative behavior of our result. First we note that for not very lar ge r apidities, 



Y ~ 1/a, we can put denominators in Eqs. ( p0|) and ( |22[ ) to be equal to 1, which through Eq. ( |12[) provides 
us with 

- dN °^ Yo) « (a P - l) 7 2 e (^-i)(^o) + S{Y Y )N*(Y ). (23) 

OYq 

The first term on the right hand side of Eq. (E3h is the usual result of the lowest order approach employing 
only one triple pomeron vertex [§,[l5| (see Fig. g). The second (delta- function) term in Eq. (^3|) corresponds 
to the contribution of the elastic cross section when the rapidity gap fills the whole rapidity interval. 

Remember that N D (Y, Yq) is the cross section of having a rapidity gap greater or equal to Yq. If one wants 
to obtain the cross section for a fixed size of the rapidity gap one has to differentiate N D with respect to Yq 
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and take a negative of the result, since N D is obviously a decreasing function of Yo. This is what was done 
in arriving at Eq. (p3|). That way we have verified that the toy model solution of our equation (j^) maps 
onto the old and well established triple pomeron vertex result Q . 

At the lowest order 7 ~ a 2 A 1 / 3 . Thus the first term on the right hand side of Eq. (p3|), which is usually 
calculated in the perturbative approaches pq| , is of the order of a 5 C 1. We note that using the nonlinear 
evolution equation ^ we were able to obtain control over the kinematic regions where the coupling is 
still small but the cross sections can become very large, with No ~ 1. Thus the traditional perturbative 
calculation of Eq. (^3|) appears like a small perturbation compared to our large cross section result of Eq. 
(ptj|), which was still obtained in the small coupling constant limit. 

At very high energies, when Y 00, one can see from Eqs. (E0h,(E2) and ([L2h that 



N D (Y,Y ) -v 1, Y -> 00. (24) 

This result is easy to understand. As we know from quantum mechanics at very high energies, when the 
total cross section is black, only the elastic and inelastic contributions to the cross section survive to give half 
of the total cross section each. The total cross section becomes a to t = 2ttR 2 , whereas totally inelastic and 
elastic cross sections will be <Ti ne i = a e i = ttR 2 . All the intermediate contributions with finite size rapidity 
gaps covering fraction of the total rapidity interval go to zero. Since N D (Y, Yo) is the cross section of having 
rapidity gap greater or equal than Yq for any finite non-zero Yo it includes the contribution of the elastic 
cross section, which corresponds to the case when the rapidity gap covers the whole rapidity interval. At 
the same time N (Y, Yo) does not include the totally inelastic contribution, when there is no rapidity gap 
at all. Thus as at Y — > 00 only the elastic piece in N D survives to give half of the total cross section, which 
corresponds to N D — > 1 in our notation, where everything has to be multiplied by ttR 2 . 

Let us define the cross section of single diffractive dissociation with the fixed rapidity gap Yq, which, for 
Yo < Y, in the light of the above discussion is 

w kvvm dN D (x 01 ,h,Y,Yo) ar(x 01 ,b,y,y ) , w vw2 , Y r9 - 

i?(xoi,b,r, Yo) = ^7 = ^7 M(Y - Y )A^ (x i,b, Y ). (25) 

OYo OYo 

In the toy model we have 

7 2 (1 - 7) 2 (qp - l) e ^p-l)(Y+Y ) 

R(Y, Y Q ) = 7 li 1) {op yje _ Y )N 2( Y y (26) 

[1 + 2 7 (e^r-VY - 1) + 7 2 (1 - 2e(^" 1 )^ + e ( aP -i)(Y+Y )^ 2 

Now one can see that the triple pomeron vertex result of Eq. (^3|) could be recovered from the first term of 
Eq. ( p6| ) by putting the denominator to be 1 and neglecting 7 with respect to 1 in the numerator. Also the 
quantum mechanical result discussed above is illustrated explicitly: for any finite Y the cross section of a 
fixed size gap R(Y, Y ) goes to zero as Y — > 00, as could be derived from Eq. ([2(f) 

Finally, let us plot the cross section of diffractive dissociation i?(Y, Yq) of Eq. (|26|) as a function of the size 
of rapidity gap Y for a fixed center of mass rapidity Y excluding the elastic cross section (delta-function) 
contribution. This is equivalent to plotting the cross section as a function of the invariant mass of the 
particles produced in the virtual photon's decay, since Yo = Y — In M\ /Q 2 (see Fig. ||) . The plot is shown 
There we put ap — 1 = .5, 7 = .03, Y = 10 and we divided Yo by Y on the horizontal axis. 
9] demonstrates that the cross section increases with the size of the rapidity gap over most of the 
rapidity interval. This is what one would expect from the lowest order formula (p3[). Also that implies that 
it is more advantageous to produce particles with smaller invariant mass M^- However at some very large 
size of the rapidity gap the cross section i?(Y, Yo) reaches a maximum and starts decreasing. The maximum 
of R(Y, Yo) in Eq. (g6) is reached at the size of the rapidity gap 

Y — = — L_ i n ( 2 SLUl + fi-^VcaP-V) . (27) 
ap - 1 V 7 7 / 

From Eq. ( jjjj ) we can immediately conclude that if the target nucleus is very large so that each of the 
dipoles undergoes many rescatterings making 7 « 1 the maximum will be pushed all the way into the 
region of negative Yo, making the cross section of Fig. [)] just a decreasing function of Yo for all rapidities 
of interest. Thus, if the effects of multiple rescatterings of individual dipoles become important then the 
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qualitative behavior of the cross section of diffractive dissociation will change completely. The effect of 
multiple rescattering is to push the maximum of Eq. (^) towards the smaller values of rapidity, which are 
easier to reach experimentally. 




FIG. 9. The cross section of diffractive dissociation in units of tvR 2 at fixed Y as a function of Y"o, which is scaled 
by Y. The elastic contribution is not included. 

If the multiple dipole rescatterings are not very important we can just use the lowest order (one rescat- 
tering) 7 ~ a 2 A 1 / 3 <C 1 and consider two limits. At not very large values of rapidity, Y ~ 1/a, when the 
single BFKL pomeron exchange is most important the maximum can never be reached since 

Y o max r ^ » — ~~r • (28) 

ap — 1 or A 1 / J ap — 1 

As energy increases and becomes of the order of saturation energy, corresponding to rapidities of the order 
of 

Y sat L^n— L- (29) 

ap — 1 a z A L i 6 

the maximum will be reached and the effect should be observed. 

We can conclude that experimental observation of the maximum of the cross section of diffractive dissoci- 
ation at certain size of the rapidity gap would signify the presence of multiple rescattering effects either in 
multiple pomeron exchanges, which are encoded in the branching of the dipole wave function, or in multiple 
interactions of each of the dipoles with the target, which were enclosed in the initial condition (|3|). Therefore 
it would be an independent test of whether the evolution equations describing the hadronic or nuclear struc- 
ture functions should be still used at the linear order or non-linear effects like in Eq. (|^) should be taken 
into account. The qualitative conclusion will also be independent of the non-perturbative initial conditions 
of the evolution equations. Here we have to warn the reader that our results in this section are based on the 
toy model and a more elaborate numerical analysis of Eq. (||) is needed to achieve the required level of rigor 
in these conclusions. Also Eq. (Q) resums fan diagrams, leaving out the pomeron loop contributions. Even 
though pomeron loops are always suppressed by powers of A compared to fan diagrams, they will become 
important at the rapidities of the order of Y ~ a _ x In ■ This is still above the saturation energy of Eq. 
( p9| ) and does not affect much the kinematic region considered above. Nevertheless for real life nuclei the 
suppression of pomeron loop diagrams is not that large and they may play an important role already at the 
rapidities of interest. 



V. CONCLUSIONS 



In this paper we have derived the small- a: evolution equation for the cross section of single diffractive 
dissociation (^). It resums multiple BFKL pomeron exchange diagrams contributing to the diffractive cross 
section. The equation is not derived in the framework of any model of multiple exchanges: it was directly 
derived from QCD. We made use of the fact that in DIS with large Q 2 the strong coupling constant is small. 
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We also employed the leading logarithmic approximation resumming logarithms of Bjorken x JIH- Finally 
we employed the large N c limit of QCD at high energies — Mueller's dipole model 

The obtained equation ^ is nonlinear and most likely can not be solved analytically. We solved a 
simplified model of the equation (Eqs. (20) and (p2[)), in which the transverse coordinate dependence was 
suppressed. We believe that this approximation preserves the qualitative features of the solution. In this toy 
model we observed an interesting phenomenon: the diffractive cross section has a maximum as a function 
of the rapidity gap Yo (Fig. ^|). This effect can not be obtained from the usual single triple pomeron vertex 
approach of [p5|p| and if experimentally observed would signify the importance of non-linear effects in QCD 
evolution of structure functions. The effect should be seen in the physically measurable quantity — DIS 
diffractive cross section and it would be an independent evidence of the onset of saturation. 
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